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have been obtained within the scope of general relativity. In particular, we 
considered Freedman-Robertson- Walker (FRW) space-time as an external ho- 
mogenous and isotropic gravitational field whereas the inhomogenous and 
isotropic Universe is given by the Godel model. The solutions obtained have 



^ ■ been thoroughly studied for different types of interaction term. It has been 

5-H 



shown that in FRW space-time equations with different interaction terms may 
have stable solutions while within the scope of Godel model only the droplet- 
like configurations may be stable, if they are located in the region where 
g 00 > 0. 
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Scalar Field with Induced Nonlinearity: 
1. INTRODUCTION 



Development of general relativity (GR) and quantum field theory (QFT) leads to the 
increasing interest to study the role of gravitational field in elementary particle physics. 
To obtain and study the properties of regular localized solutions to the nonlinear classical 
field equations is motivated mainly by a hope to create a consistent, divergence-free theory. 
These solutions, as was remarked by Rajaraman TO] give us one of the ways of modeling ele- 
mentary particles as extended objects with complicated spatial structure. In such attempts 
it is natural to treat the field nonlinearity not only as a tool for avoiding the theoretical 
difficulties (such as singularities) but also as a reflection of real properties of physical sys- 
tem. It should be also emphasized that the complete description of elementary particles 
with all their physical characteristics (e.g., magnetic momentum) can be given only in the 
framework of interacting field theory [Q . In this paper we present some regular particle-like 
solutions within the scope of general relativity for an interacting system of scalar and elec- 
tromagnetic fields, confining ourselves to static, spherically and/or cylindrically symmetric 
configurations. 



2. FUNDAMENTAL EQUATIONS 

We consider a system with the Lagrangian 

where the first term describes the Einstein gravitational field (R is the scalar curvature, 
k = 8irG is the gravitational constant) and *ff((p) is some arbitrary function characterizing 
interaction between the scalar {ip) and electromagnetic (F^) fields. This kind of interaction 
has been thoroughly discussed in ||. The nongravitational part of the Lagrangian( [2.1|) 
describes a system with a positive-definite energy if and only if ^ > 0. Let us write the 
scalar and the electromagnetic field equations corresponding to the Lagrangian( |2.1| ) 
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'—g dx L 



dx^J 



16/T 



0, * w = — 



dip 



(2.2) 



1 d_ 

/—g dx h 



-gF v ^)) =0 



(2.3) 



The corresponding energy-momentum tensor reads 



16n 



(2.4) 



3. SPHERICALLY SYMMETRIC CONFIGURATIONS 

A. Solutions in FRW Universe 

As an external homogenous and isotropic gravitational field we choose the FRW space- 
time. This Universe is very important as the corresponding cosmological models coincides 
with observation. The interval in the FRW Universe in general takes the form 

dr 2 



ds 2 = dt 2 - R 2 (t) 



1 — kr 2 



+ r 2 {d$ z + sin 2 $d(j) 2 } 



(3.1) 



Here R(t) defines the size of the Universe, and k takes the values and ±1. We consider 
the simple most case putting R(t) = R = constant, which corresponds to the static FRW 
Universe. In static case k = corresponds to usual Minkowski space, k = +1 describes 
the close Einstein Universe and k = — 1 corresponds to the space-time with hyperbolic 
spatial cross-section. Note that the velocity of light c has been taken to be unity. 

As was mentioned earlier, we seek the static, spherically symmetric solutions to the 
equations ( |2.2|) and ( |2l|). To this end we assume that the scalar field is the function of 



r only, i.e. tp = ip(r) and the electromagnetic field possesses only one component Fio 
dA /dr = A'. 

Under the assumption made above, the solution to the equation ( [2.3|) reads 



^ =^) JL l^ir- (3-2) 
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where q is the constant of integration and P(ip) = l/^((p). Putting ( |3.2|) into ( |2.2|) for the 
scalar field we obtain the equation with "induced nonlinearity" M 



:i-*rV + L -?* ! V- ] grf, = 0, « 2 = £ (3.3) 



This equation can be written in the form 

d 2 ip 2q 2 



~ -kP* = (3.4) 



with the first integral 



dz 2 R 2 



dip 2q 



where we substitute z = y 1/r 2 — k. Here Cq is the constant of integration, which under the 
regularity condition of Tq at the center turns to be trivial, i.e., Cq = 0. Finally we write the 
solution to the scalar field equation in quadrature 

In accordance with ( pT2| ) and ( |3.5| ) from ( |2.4j) we find the density of field energy of the 
system 

TS - ^ (3.7) 

and total energy of the material field system 

E f = j T °y^rf 3 x = -8vrg / VPdip (3.8) 

Thus, we see that the energy density Tq and total energy Ef of the configurations obtained 
do not depend on the conventional values of the parameter k. As one sees, to write the 
scalar (ip) and vector (A) functions as well as the energy density (Tq) and energy of the 
material fields (Ef) explicitly, one has to give P((p) in explicit form. Here we will give a 
detailed analysis for some concrete forms of P(p)- Let us choose P(ip) in the form 

P{<p) = Pocos 2 (^) (3.9) 
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with A being the interaction parameter. Inserting ( |3.9|) into ( |3.4| ) we get the sin-Gordon 
type equation || 

d 2 <p Xq 2 P . 



dz 2 + AR 2 



-sin(Ay) = 



(3.10) 



with the solution 



(p(z) = — arcsintanh[&(,2 + zi)], b = ^ ° , Z\ = const 
A R 



(3.11) 



Let us analyze the solution ( |3.11|) . It can be shown that lim<£> = tt/X for all k. As one sees, 

r^O 



the solution possesses meaning only in the region where z = yl/r 2 — k > 0. It means, in 
case of k = +1 the configuration confines in the interval < r < 1. Then the asymptotic 
behavior of the solution ( |3.11 ) can be written as 



V 3 



0, k = —1, z\ = —1, r — > oo 
0, k — 0, z\ — 0, r — > oo 
0, fc = +l, «i = 0, r>l 



(3.12) 



From ( |3.8| ) we find the total energy of the system Ef = — 167iqy/P^/ X. For the choice of 
P((p) in the form 



2\2 



P = \(a 2 - <p 2 ) 



(3.13) 



with A being the coupling constant and a being some arbitrary constant, from ( |3~D we 
obtain 



4AaV - 4Ay? 3 + 



fe 2 







(3.14) 



The equation ( |3.14| ) can be seen as an MKdV one. Indeed, a KdV equation 



du „du „d 3 u 

— + au p — + (3— = 0, 

at ax ox 6 



(3.15) 



can always be converted to 



u 



P+i 



A 2 



a 



Du + a + p— = 

p+1 dz z 



(3.16) 
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if one looks for stationary solution of the form u = u(z) where z = x — Dt. In our particular 
case p = 2 and the equation (|3.14|) is an MKdV one. The scalar field function in this case 
has the form 

<p(z) = atanhfv Xab(z + z 2 )], b = — (3-17) 

R 



Taking into account that z = Jl/r 2 — k one sees that at the origin limy? = a whereas at 

v 1 — >o 

the asymptotic region for different value of k we get 



0, k = —1, z 2 = —1, r — > oo 

0, k = 0, z 2 = 0, r^oo (3.18) 
0, k = +1, z 2 = 0, r > 1 

From ( |3.8|) we find the total energy of the system to be Ef = — lQnq^/Xa 3 /3. A specific 
type of solution to the nonlinear field equations in flat space-time was obtained in a series 
of interesting articles 0. These solutions are known as droplet-like solutions or simply 
droplets. The distinguishing property of these solutions is the availability of some sharp 
boundary defining the space domain in which the material field happens to be located, i.e., 
the field is zero beyond this area. It was found that the solutions mentioned exist in field 
theory with specific interactions that can be considered as an effective one, generated by 
initial interactions of unknown origin. In contrast to the widely known soliton-like solutions, 
with field functions and energy density asymptotically tending to zero at spatial infinity, the 
solutions in question vanish at a finite distance from the center of the system (in the case 
of spherical symmetry) or from the axis (in the case of cylindrical symmetry). Thus, there 
exists a sphere or cylinder with critical radius tq outside of which the fields disappear. 
Therefore the field configurations have a droplet-like structure [|7|,|9|,|T0H . 

To obtain the droplet-like configuration we choose a very specific type of interaction 
function P((p) which has the form fllT 



P(<p) = J 2 " 4 /^l- J 2 '^ 2 (3.19) 
where J = \<p, a — 2n + 1, n — 1, 2 • - -. Putting ( |3.19|) into ( ^6|) for cp one gets 



v>W = - x 
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l-exp[-—(z-Zo, 



a/2 

(3.20) 



Recalling that z = Jl/r 2 - k from ( |3.20| ) we see that at r — > the scalar field tp takes 



the value <p(0) — > 1/A and at r — > r c = 1/Jzq + fc, the scalar field function becomes trivial, 
i.e., <p(r c ) — > 0. It is obvious that for r > r c the value of the square bracket turns out to be 
negative and <p(r) becomes imaginary, since a is an odd number. Since we are interested in 
real ip only, without loss of generality we may assume the value of ip to be zero for r > r c , 
the matching at r = r c (i.e., z = zq) being smooth. Note that, for k = +1, the scalar field 
is confined in the region < r < 1, as it was in previous two cases. The total energy of the 
droplet we obtain from ( |3.8|) has the form 

E > = W=T) < 3 ' 21 > 

As is seen from ( |3.21| ), the value of the total energy does not depend on the size of the 
droplet, it means droplets of different linear size share the same total energy. 

B. Stability problem 

To study the stability of the configurations obtained we write the linearized equations 
for the radial perturbations of scalar field assuming that 

<p(r,t) = <p(r) + £(r,t), (3.22) 

Putting ( |3.22j ) into ( |2~2|) in view of ( |3.3|) we get the equation for £(r, t) 
•■ R. 1 - kr 2 „ 2 - 3kr 2 . q 2 P w 

£+3 s«~s-^~^- f+ ^H t{ = (3 ' 23) 



The second term in ( |3.23| ) is zero since we assume the FRW space-time to be static one 
putting R = constant. Assuming that 

£(r,t)mv(r)exp(-iSlt), Vl = uj/R (3.24) 

from ( |3.23|) we obtain 
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2-3£t 2 r 2 q 2 P v 



(1 _ kr 2 )v" - v 1 + u 2 - v = (3.25) 

The substitution 



v(0 = r ' v ( r )i C = — ^arcsin(v / A;r) (3.26) 



leads the equation ( |3.25|) to the Liouville one |l2j 



% + (^-l%», = 0, ^) = -* + ^(_^)' (3.27) 

For the interaction term ( |3l]) we see that 

2./ \ &C 4 i2 2 /sin( v / A?C)\ 4 1 



V^)>0, if and only if tanh 2 fr(,z + z x ) > ( ^ j +^ (3.28) 

Thus we find that the equations with trigonometric nonlinearity contain stable solutions. 
Given P(ip) in the form ( |3.13| ) we find that 



* i „ , 9l . , k( 4 R 2 /sinfv^Ox 4 a 2 

V(¥>)>0, if tanh 2 fe(z + ^)> T |- iX ( ^1 ^ ) + y (3.29) 



As is seen from ( |3.29|) the equations with polynomial type of nonlinearity too contain some 



stable solutions. For the droplet-like configurations, i.e., for the interacting term P((p) given 
by ( |3.19|) , it can be shown that the potential 



liml%) -»■ +oo, Jim V(<p) -> +oo (3.30) 

beginning with a > 5. It means that the droplet-like configurations ( |3.20| ) with a > 5 are 
stable for the class of perturbation, vanishing at r = and r = r c . 

4. CYLINDRICALLY SYMMETRIC CONFIGURATIONS 
A. Solutions in Godel Universe 



In the previous section we studied the possibility of formation of regular localized con- 
figuration in homogenous and isotropic FRW Universe. Let us now continue our study in 
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the homogenous but anisotropic Universe. In particular we consider the model proposed by 
Godel. The linear element of Godel Universe in cylindrical coordinates reads III 



ds 2 = dt 2 - dp 2 + -Usinh 4 fip - smh 2 n P }d4> 2 - ^smh 2 Qpd4>dt - dz 2 (4.1) 

where the constant Q is related with the angular velocity u: u = \/2Vt. This form of 
linear element of the four- dimensional homogenous space S directly exhibits its rotational 
symmetry, since the g^ u do not depend on <fi. It is easy to find 

lim J — q = lim — -sinh(2f2p) — ► p 

i.e. at uj —>■ Godel Universe transfers to the flat one. 

As in spherically symmetric case, here too we seek the static solutions to the equations 
( |2.2|) and ( |2.3|) assuming the scalar field to be the function of p only, i.e., ip = <p{p) and the 
electromagnetic possesses only one component Fio = dA /dp = A'. For the electromagnetic 
field in this case we find 



2QDP/smh(2Qp), D = const (4.2) 



The scalar field equation ( |2.2| ) with regards to ( |4.2[ ) reads 



Putting y = ^lntanh(fip) from ( |4~3| ) one gets 



dy 2 



8q 2 n 2 P^ = (4.4) 



with the first integral 



^ = ±4^0^ + D (4.5) 

Here Do is the constant of integration, which under the regularity condition of Tq at the 
center turns to be trivial, i.e., D = 0. Finally we write the solution to the scalar field 
equation in quadrature 
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dp 



p 4gfi(y - y ) (4.6) 
In accordance with ( |4.2|) and ( |4.5| ) from ( |2.4j ) we find the density of field energy and the 



total energy of the system 



_ 10qWP_ 

J- n — _ o , 



sinh 2 (2fip) (4 - 7) 
E f = 8nq J VPdp (4.8) 

As in the previous case, to write the scalar (p) and vector (A) functions as well as the 
energy density (Tq ) and energy of the material fields (Ef) explicitly, one has to give P(p>) 
in explicit form. Here again we will thoroughly study the solutions obtained for different 
concrete forms of P(p)- Choosing P(<p) in the form ( |3.9|) , i.e., 

P(p) = P cos 2 (^) (4.9) 



with A being the interaction parameter, from ( O ) we get the sin-Gordon type equation 



d 2 cp 



dy 2 

The solution to this equation can be written in the form 



+ 4Xq 2 tt 2 P sm(Xp>) = (4.10) 



<P(P) = \ 



arctan 



f tanhQp\ a n 



\ flpo ) 4_ 

where po is the constant of integration, giving the size of the system. Without losing the 
generality we can choose a > 0. Then one finds lim<p = —tt/\. For p > the field 
p steadily increases up to tt/X. In particular, at spatial infinity we get lim ip = 0. In 
this case P(poo) = 1 which corresponds to the exclusion of interaction at spatial infinity. 
The total energy of the system in this case coincides with that of in FRW Universe, i.e., 
Ef = — 16Tcq^/P ) / \. Choosing P((p) in the form ( |3.13| ), i.e., 



a = XqJP (4.11) 



P = X(a 2 - p 2 ) 2 (4.12) 
from ( |4.4|) we as in FRW case again obtain MKdV type equation with the solution 
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tanh a (fip) - 1 r- . 

From ( [4.13Q it is clear that lim<p — > —a and lim ip —>■ From ( |4.8| ) we find the total 
energy of the system to be Ef = — lQnqy/Xa 3 /3 as it was in FRW case. The choice of the 
interaction term P((p) in the form ( |3.19| ), i.e., 

P{ip) = J 2 - A /°(l - J 2/fT ) 2 (4.14) 

where J = \<p, a — 2n + 1, n — 1, 2 • • •, leads to the following expression for scalar field 

/ tanh(fip) 



<p(p) = \ 



, a = (4.15) 
a 



^tanh(fipo 

where po is an arbitrary constant. For a > the solution possesses physical meaning at 
p < po and becomes meaningless at p > po- Outside the cylinder p = po one can put <p = 0. 
This trivial solution is stitched with the solution at p = p under condition <p'(po) = 0, 
which fulfills if and only if 4A|g| > a. Consequently, at p > p the Lagrangian becomes 
physically meaningless, however its limiting value at p — ► po — is equal to zero. Continuing 
it at p > po, one can consider the field be totally trivial in this area. Thus we get the 
droplet-like configuration. The field <p steadily decreases from <p(0) = 1/A to <p(po) = with 
<p'(0) = (for a > 3) and <p'(Po) = 0. The total energy of the "droplet" is defined as 

which remains unaltered even in flat space-time (Q = 0). It means that the "droplet" does 
not feel Godel gravitational field. 



B. Stability problem 



Let us now study the stability of the configuration obtained. In doing so we consider 
the perturbed scalar field 5<p = x(p, t) that leaves the cylindrical-symmetry of the system 
unbroken. The linearized equation for the perturbed scalar field looks 
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' (j / ^grxr) + J—.P^x = (4.17) 



v^fe" v ' ; 2{-g) 

Since, for the case considered, g 00 = (1 — sinh 2 f2p)/cosh 2 f2p, it is clear that the type of 

equation (|4.17f) changes on the surface where sinhfip = 1. The Cauchi problem for ( f4.17|) 

is incorrect by Hadamard in the region, where g 00 < 0. In connection with this only the 

droplet-like solutions can be stable by Lyapunov, if they are located in the region where 

g 00 > @ . Assuming that 

X(p,t) = v(p)exp(-iet) (4.18) 

from ( |4.17| ) we get 

_ , . „ . , r l-sinh 2 fip 9 2Vt 2 q 2 

The substitution 

T}(0 = Z(p) ■ v(p) (4-20) 

where 



c = / v \Jhfy" P ( fo £ = t 4sinh2 - sinh2 ^)] 1/4 

leads the equation for perturbed field to the Liouville one 

^ + [z 2 -V]v = (4.21) 

Here the effective potential V(ip) takes the form 

tt 2 r q 2 P w (4sinh 6 Qp- 16sinh 4 fip + 3sinh 2 fip- l)cosh 2 Op. 
^ ~~ (1 - sinh 2 fip) ^ 2sinh 2 fip + 4sinh 2 fip (1 - sinh 2 fip) ' 

For the interaction function, chosen in the form ( |3.19| ), one finds 

x2 > 2 -12a + 16 4a 2 -12a + 8 



(7 



Taking into account that lim J — > 1 and lim J — > 0, it can be shown that 

p— >0 p— >po 

lxmV((p) — > +00 for cr < 4gA 

p^O 



lim V(<p) — ► +oo for a > 5 

p->po 

Here we used the fact that g 00 > 0, i.e., sinhfip < |1|. Thus, as in the previous case we find 
that the droplet like configurations are stable for 5 < a < 4|g|A. 



B. Saha 



13 



5. CONCLUSION 

We obtained the regular particle-like solutions to the scalar field equations with induced 
nonlinearity in external gravitational fields described by Freedman-Robertson- Walker and 
Godel Universes respectively. Beside the usual solitons, a special type of regular localized 
configurations, known as droplets, have been obtained. It has been shown that the droplet- 
like configurations possess limited energy density and finite total energy and the droplets of 
different linear sizes up to the soliton share one and the same total energy. It is noteworthy 
to notice that in the spherically symmetrical case (i.e., in the FRW Universe) at r c — > oo 
for k = droplet transfers to usual solitonian solution, while for k — ±1 this is not 
the case. It has also been shown that in FRW space-time equations with different type 
of nonlinearities may contain stable solutions, whereas in case of Godel Universe only the 
droplet-like configurations may be stable. It is noteworthy to remark that in FRW space- 
time with k — +1 the field function is confined in the region < r < 1 independent to the 
choice of interaction function P(<p). 
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FIGURES 



FIG. 1. Perspective view of droplet-like configurations for different values of a. Here the 
thick-line, dash-line, dash-dot-line and thin-line correspond to the value of a = 3, 5, 7, 9 respectively. 
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FIG. 2. Perspective view of interaction function P((p) providing droplet-like configurations 
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